Long-range resonant interactions in biological systems 
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The issue of retarded long-range resonant interactions between two molecules with oscillating 
dipole moments is reinvestigated within the framework of classical electrodynamics. By taking 
advantage of a theorem in complex analysis, we present a simple method to work out the frequencies 
of the normal modes, which are then used to estimate the interaction potential. The main results 
thus found are in perfect agreement with several outcomes obtained from quantum computations. 
Moreover, when applied to a biophysical context, our findings shed new light on Frohlich's theory 
of selective long-range interactions between biomolecules. In particular, at variance with a long- 
standing belief, we show that sizable resonant long-range interactions may exist only if the interacting 
system is out of thermal equilibrium. 
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Introduction — Living organisms are well-known to sup- 
port many biochemical processes which, besides being ex- 
tremely specific, seem to follow a precise time schedule. 
In this context, the motion of the molecules can hardly 
be described on the basis of thermal fluctuations only. 
Most studied concepts to this effect refer to electromag- 
netic short-range interactions [l[ , which leaves long-range 
interactions poorly investigated thus far. This is essen- 
tially because in biological systems free ions of cell water 
tend to screen any electrostatic potential at a distance 
that usually does not exceed a few angstroms. However, 
this screening proves generally inefficient for interactions 
involving oscillating electric fields. To this respect, in 
1972, Frohlich emphasized Q that two molecular sys- 
tems which would exhibit large oscillating dipole mo- 
ments could interact via a long-range potential, provided 
that the oscillation frequencies were roughly similar. 

Applied to biological systems, it was suggested that 
such resonant forces would have a profound influence on 
the displacement of specific biomolecular entities, and 
thus on the initiation of a particular cascade of chemical 
events. Later on, possible examples of such interactions 
were reported - at the cellular level - between erythro- 
cytes [3j. In parallel, it was shown that the membranes 
of these cells have the ability to oscillate at frequencies 
of about 40 GHz jij. Moreover, since then [f|, numer- 
ous evidences of electromagnetic long-range interactions 
between cells have been identified including yeast cells 
Very recently, supported by experimental evidences 
of collective oscillations in biomacromolecules (see Refs. 
in 0|), possible experimental tests were proposed to as- 
sess whether such interactions could be relevant at the 
biomolecular level 0|. 

Coming to Frohlich's predictions, an intriguing result 
of his computations is the possibility of observing long- 
range resonant interactions even if the system of oscillat- 
ing dipoles is at thermal equilibrium. This would occur 
when the retardation of the electric field mediating the 



interaction is relevant. For large intermolecular distances 
r, the potential was given at thermal equilibrium as 
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where e(w±) is the permittivity of the medium at the 
normal frequencies tu± of the interacting system. 

Though Frohlich's theory has been well received and 
developed for a long time by many authors Q , a thorough 
derivation from first principles of his main results on res- 
onant interactions, as for instance the potential given in 
Eq. |T]), is still lacking. In anticipation of advanced ex- 
perimental investigations 0| , the present paper is thus in- 
tended to be a more comprehensive and analytic account 
of this theory and to put the mentioned interactions at 
their right place within the framework of electrodynamic 
forces. To that purpose, we calculate the normal frequen- 
cies of a system of two interacting oscillating dipoles in 
an exact way, including field retardation. This is made 
possible by means of a theorem in complex analysis, the 
Lagrange inversion theorem, that is detailed below. We 
then show how these results can be used to estimate the 
interaction energy between two atoms on one hand, and 
two oscillating dipoles on the other hand. In the former 
case, a connection with known QED results is provided 
and in the latter one the connection with Frohlich's the- 
ory is made. In particular, we prove that the form for the 
potential |lj is misleading, suggesting that long-range in- 
teractions exist only if the system is out of equilibrium. 

Theory — Let us consider two molecules A and B with 
oscillating dipole moments fx A and fi B . The equations 
of motion can be given in general terms as 

j fj-A + lAfJ-A+^A^A = (AE B (r A ,t) + f A {H A ,t) 

I Ms +7sMs + uj 2 B fi B = ( B E A (rB,t) + f B (fi B ,t), 

(2) 

where uja,b and ja, b are the harmonic frequencies and 
damping coefficients of the dipoles. Here, the interaction 



2 



takes place through the electric field Ea,b(t, t) generated 
by each molecule, located at r = r B ,A, while the coupling 
constants are given by £a = Q A /mA where Qa and tua 
are the effective charge and mass involved in the dipole 
A, and similarly for Q B . Finally, f A and f B are func- 
tions accounting for possible anharmonic contributions, 
thermal noise as well as possible external excitations. 

Normal modes analysis — A possible way of estimating 
the interaction energy of the system (carried out in the 
last part of this letter) is to use the normal frequencies 
of the associated harmonic conservative system given by 



ft A + UaPa = (AE B (r A ,t) 



sMb = CBE A (r B ,t). 



(3) 



The normal frequencies are defined as the frequencies 
u) such that /j, a b(*) = t^A B e are solutions of ((3|). In 
this context, the electric field generated by a harmonic 
dipole is well-known from classical electrodynamics fill ) 
and one has the following matrix relation 



E B (r A ,t) = x{r,w)ii B t 



(4) 



with a similar expression for EA(r B ,t). Here, x re P~ 
resents the generalized susceptibility matrix (retarded 
Green function) of the electric field - note that the spa- 
tial dependence of x hi ([!]) reduces simply to the distance 
r = \ta — Tb\ because the intermediate medium is sup- 
posed to be homogeneous and isotropic. In addition, to 
model an aqueous biological environment filled with free 
moving ions, in what follows we will consider it as a dis- 
persive dielectric medium of complex permittivity s(uS). 

At this stage, it should be remarked that the general- 
ized susceptibility of the electric field is analytic every- 
where on the complex plane except on an uncountable 
subset where it has a discontinuity. In a homogeneous 
dielectric, this subset is given by {ui <E C | lm(oj^/e(uj)) = 
0} while x reads as [ll[ 



Xii(r,cj) = X22(r,uj) 



e(uj)i 



X33[r,u) = -— 1 =F — 

£(w)r c> I c 
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(5b) 



and Xij{ r i w) = when i ^ j\ the ± sign is attributed to 
positive or negative values of Im(w^/e(w)), respectively. 
Let us also specify that the diagonal form of x is due to 
the choice to set the z axis along r. Finally, it should be 
stressed that for real values of ui, each \a is a complex 
number whose imaginary part essentially accounts for the 
dissipation due to the field propagation jl2j |. Bearing in 



mind that in computing normal frequencies, one drops 
dissipation effects, only the real parts denoted by x'u will 
be thus considered in what follows. 

Substituting into Eq. ([3]) the expected harmonic forms 
for fi A b and Eq. ^ , one obtains a system of linear 
algebraic equations. The existence of the solutions is 
then ensured by the vanishing of the determinant (w^ — 
cj 2 )(o;| - lj 2 ) - CaCb(Xi,(^ w)) 2 = 0. This indicates two 
possible values for u> > for each i, namely, u>i,+ and 
verifying 




= g V % J + ^ B( - X '^ r,UJi ' ± ^ 2 - 

(6) 

At this point, it must be emphasized that each normal 
frequency is strongly dependent on the proximity of the 
frequencies of the dipoles, as detailed below. 

Resonant case — When lo a — = <-^B, equation ([6]) 
is readily simplified and at large separations, one can 
expand u)i± around cjq as 



with 
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The implicitness of the equation can then be solved 
by making use of the Lagrange inversion theorem [l3| : 
let be a contour in the complex plane surrounding a 
point a, and let <f> a function analytic inside and on c £ '. 
If t € C is such that the inequality \t(f>(z)\ < \z — a\ is 
satisfied for all z on the perimeter oftf, then the equation 
lj = a + t(j>(u>) in ui has one root inside and any further 
function g analytic inside and on can be expanded as 
a power series in t by the formula 



g(u) = g(a) + 



t n d n 



^ n\ da n 



n=l 



dg 
da 
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Identifying t<j){uj) with <fi^±(r,u>) of equation ([7]), one 
must carefully choose the contour ^ out of the domain 
of discontinuity of x'u ■ On the other hand, since each x'u 
(and so </>,) is a sum of inverse power laws of r, we can 
assume that for large enough separation the inequality 
4>i.±{ r i z ) < \ z ~ w o| holds for z close enough to loq so 
that 'if may be still chosen in the domain of analyticity 
of x'u- Letting g as the identity function and n < 2, one 
finds after some algebra 



w«,±(r) ~ loq ± 



VCAtBX'u(r, ujq) 
2w 

Ca(b d 



2 du 



(9) 



to second order in x'u 



1,2,3. The normal frequen- 



cies thus found are equal to the resonance frequency uiq 
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plus a shift due to the (effective) interaction. The first 
contribution of each frequency shift is proportional to 
the real part of the susceptibility matrix elements. As 
detailed below, this term is responsible of the long-range 
nature of the resonant interaction energy as, according to 
([5]), each x-element reads as a polynomial in l/r Q with 
a < 3 (the dimension of physical space) : in the limit 
r <C c/luo (near zone limit), it goes as ±l/r 3 with the 
intermolecular distance while it oscillates at longer dis- 
tance (intermediate and far zone limits) with a ±1/V 2 or 
even ±l/r envelope. 

Off-resonance case — On the contrary, when uia 3> Wb 
(or similarly when uia *C ub), one can use Eq. (j6]) 
to approximate u>i,± around usa.b as u>i,± ~ uia,b ± 
GiCs (Xri(f,w ij± )) 2 / \2wa,b{ua - u>%)] to second order 
in Xu- Again, for large r, the Lagrange inversion theorem 
can be applied so that tip of the theorem corresponds to 
the last term of the equation. To compare the results 
with those of the resonant case, it is enough to apply Eq. 
([5]) up to n = 1. Thus, the normal frequencies are 



,±(r) ~ wa,b ± 



CaCb (Xu(r, ua,b)Y 
2uj a ,b (u 2 A - u%) 



(10) 



At variance with the resonant case, the frequency shifts 
associated with the unperturbed frequency uja.b are now 
proportional to (x'ii) 2 a t first order. When r <C c/uja.b, 
this leads to a short-range ±l/r 6 contribution while at 
very long distances it oscillates decaying as ±l/r 2 . 

Interaction energy ( atoms ) — Now that normal frequen- 
cies have been worked out, it is interesting to remark 
that the energy shifts U±(r) = hj^i ( w i,±( r ) — wa,b) are 
identical with the real photon contributions that appears 
in the interaction energy of two two-level atoms in an ex- 
cited state, as it can be found within the framework of 
QED. Since the atoms and the radiation field mediat- 
ing the interaction can be considered as an ensemble of 
harmonic coupled oscillators, the correspondence is not 
so surprising if we add to this that normal modes and 
susceptibilities are the same for classical and quantum 
oscillators. Thus, when uja — <-^o = ^b, U± computed 
from Eq. ((9|) are equal to the energy shifts due to the in- 
teraction between two atoms with a common transition 
frequency ujq [HI- Alternatively, in an off-resonance sit- 
uation, one finds identical expressions for U± computed 
from Eq. (JTUJ) and the energy shifts of two atoms with 
distinct transition frequencies, uja and ujb respectively 
(with loa > wb here) [15j . Readers interested in a more 
detailed investigation through the references mentioned 
above will identify (a as um|/M.i| 2 for all i = 1,2,3, 
where jiA.i is the transition dipole moment of atom A 
along the i spatial coordinate (the same for atom B) . 
The polarizability of each atom is then simply given by 



QA,s(w) 



A,B 



As a complement to these results, two important points 
should be mentioned. First, the above classical compu- 
tation of normal frequencies does not allow to deduce 
the energy contribution due to virtual photons that ac- 
counts for the interaction between the ground states of 
the atoms. This is confirmed by the fact that this con- 
tribution which is always present whether the system is 
excited or not, was shown to arise purely from vacuum 
fluctuations Il6|. Secondly, if we refer to recent quantum 
calculations |17| . the real photon contribution of the in- 
teraction energy in the case uja ^b would be noticeably 
different from the one suggested in the present paper as 
well as in the references mentioned in [15(. This result is 
distinguished by the presence of \xu\ 2 instead of (x'u) 2 in 
Eq. (fT0|) . that leads to a spatially monotonic - instead of 
oscillating - potential. Nevertheless, deeeper theoretical 
investigations [l8| revealed that the monotonic potential 
holds true on time scales much larger than the sponta- 
neous decay time of the excited atom(s). In the opposite 
case, the oscillating potential is obtained, thus recovering 
the above mentioned quantum-classical correspondence 

Interaction energy (real dipoles) — Contrarily to the 
case of two atoms, the interaction between real oscillating 
dipoles, i.e., molecules whose oscillating dipole moments 
are not due to electron motions but rather to confor- 
mational oscillations, have not been given much atten- 
tion in the literature. As mentioned in the Introduction, 
a remarkable exception is given by Frohlich who tack- 
led this issue forty years ago in a biophysical context. 
At that time, he emphasized in particular that resonant 
long-range interactions may occur between two harmonic 
dipoles even if none of the normal modes is excited be- 
yond thermal equilibrium. However, we show that this 
statement is actually incorrect. 

To clarify, let us consider the system of dipoles A and 
B in thermal equilibrium and suppose as an example that 
fnoA,B "C ksT so that classical effects are dominant [20| . 
The interaction energy is then given by the difference of 
free energy of the coupled and uncoupled systems, i.e., 
U(r) = F(r) - F(oo) = -k B Th\[Z{r)/Z{oo)}, where 
Z(r) is the partition function of the system when the 
dipoles are separated by a distance r. Using Boltzmann 
distribution for each normal mode, one can easily show 
that the interaction energy has the following form Q 



U(r) = k B T^2\n 



LOALOB 



(11) 



Taking loa — = w b> the result exposed by Frohlich 
is obtained by substituting w,.±(r) and u)i t —(r) for their 
implicit form (|T|). Then, at first order, U reads as 



U(r) 



J A,B 
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In addition, Frohlich only considered the limit r <C 
c/ljq together with the condition e ~ e' = Re(e) (non 
absorbing medium in the considered frequency range). 
From Eq. ([5]), one has in this case : 



k B Ty/( A (B 1 \ - 

u (r)~ — —2 — isZ^ 



1 



1 



£'(wj,+) e'(wi-) 



(12) 

with (7i=c*2 = — 1 and 173 = 2. At this stage, Frohlich 
claimed that the above 1 /r 3 form for U may account for 
the existence of long-range resonant interactions in ther- 
mal equilibrium provided e (k>i,+) 7^ £ However, 
as already mentioned, this form arises simply since the 
implicitness has been not resolved yet. Hence, by using 
the Lagrange inversion theorem and by substituting, for 
example, g of equation ((SJ for 1/e', one get immediately 
that the terms in curly brackets in Eq. (|12[) vanish at 
first order. As a result, expansion to second order shows 
that U is actually proportional to 1/r 6 . To compute the 
complete form for U, one then needs to come back to Eq. 
([lip . Using the explicit form of the normal frequencies 
derived above ((9]), one easily obtains 



U(r) 



3k B T( A ( B 1 



d\n[e'(oJo)] 
duo 



(13) 



where the ability of the potential to be attractive or 
repulsive depends on the derivative of e at uiq, represen- 
tative of the dispersive properties of the medium. 

In the end, it goes without saying that the long-range 
contributions that appear in the expression of the normal 
frequencies at resonance, Eq. ©, cancel each other since 
the energy of both modes is given by Boltzmann distri- 
bution; this happens independently of considering retar- 
dation effects. By inference, long-range resonant interac- 
tions between two dipoles (biomolecules) may occur only 
beyond thermal equilibrium provided that the excitation 
of one normal mode is statistically "favored" compared 
to the other. In this case, anharmonicity of the dipoles, 
as mentioned in commenting Eqs. would allow en- 
ergy redistribution among normal modes whereas energy 
supply could be essential to maintain a high degree of 
excitation despite energy losses. In a biological context 
this energy supply may be attributed to environmental 
metabolic activity. This scenario was depicted by Froh- 
lich in general terms. He showed that a set of coupled 
normal modes can undergo a condensation phenomenon 
characterized by the emerging of the mode of lowest fre- 
quency containing, in the average, nearly all the energy 
supply [Ioj|- Here such a process would result in an ef- 
fective attractive potential whose amplitude is dependent 
on the "stored" energy Q (of course, the above given re- 
marks still apply in the quantum case when Hujq 3> k B T). 
To conclude, it is worth noting that retardations effects 
at large r bring about interactions with a 1/r dependence 
[last terms in Eqs. (|5|)]. i.e., of much longer range with 



respect to the interactions proposed by Frohlich. This 
last result could be of utmost relevance for a deeper un- 
derstanding of the highly organized molecular machinery 
in living matter, as emphasized in the Introduction. 
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